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Abstract 

We analyze the role played by the Polyakov loop in the description of the chiral phase transition 
within the framework of nonlocal SU(3) chiral models with flavor mixing. We show that its presence 
provides a substantial enhancement of the predicted critical temperature, bringing it to a better 
agreement with the most recent results of lattice calculations. 
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The detailed knowledge of the phase diagram for strongly interacting matter has become 
an issue of great interest in recent years, both from the theoretical and experimental points of 
view. On the theoretical side, even if a significantprogress has been made on the development 
of ab initio calculations such as lattice QCD , these are not yet able to provide a full 

understanding of the QCD phase diagram due to the well-known difficulties of dealing with 
finite chemical potentials. In this situation, it is important to develop effective models that 
show consistency with lattice results and can be extrapolated into regions not accessible by 
lattice calculation techniques. In previous works 

mm 

the study of the phase diagram of 
SU(2) chiral quark models that include nonlocal interactions js] has been undertaken. These 
theories can be viewed as nonlocal extensions of the widely studied Nambu— Jona-Lasinio 

CI 

model [9|. In fact, nonlocality arises naturally in the context of several successfu l ap proaches 
to low-energy quark dynamics as, for example, the instanton liquid model lOj and the 
Schwinger-Dyson resummation techniques 11 1. Lattice QCD calculations |l2| also indicate 
that quark interactions should act over a certain range in momentum space. Moreover, 



several studies 



13 



141 ] have shown that nonlocal chiral quark models provide a satisfactory 



description of hadron properties at zero temperature and density. On the other hand, 
when looking at the description of the chiral phase transition, it has been noticed that for 
zero chemical potential these SU(2) models lead to a rather low critical temperature Tj?^ in 
comparison with lattice results [4,5]. The aim of the present work is to go one step beyond 
these previous analyses, studying the finite temperature behavior of nonlocal chiral models 
that include mixing with active strangeness degrees of freedom, and taking care of the effect 
of gauge interactions by coupling the quarks with the Polyakov loop. The inclusion of the 



Polyakov loop has 



PNJL models) 
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3een considered recently in the context of NJL-like models (so-called 
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17 
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19|, serving as an order parameter for the deconfinement 



transition . In particular, this has been done in Ref. 20] in the framework of a nonlocal 
two-fiavor model, focusing on the analysis of mesonic correlations. Here we show that within 
nonlocal SU(3) effective models with flavor mixing the coupling to the Polyakov loop leads 
to a significant increase of the chiral phase transition temperature T°^, which is otherwise 
as low as in SU(2) symmetric models. As stated, this is a desired effect in order to be in 
better agreement with lattice expectations 211 ]. 

We deal here with the nonlocal covariant SU(3) quark model described in Ref. 22| . 
including now the coupling to the Polyakov loop. The Euclidean effective action for the 
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quark sector of this model is given by 

Se^ J d'x [-zj,d, + m] i^ix) - I [j!{x) j!{x)+j:ix) j^ix)] 

-f T^bc [fa{x)j^{x)j!{x) - 3 j!{x)jax)j^{x)] } , (1) 

where ip is a chiral U{3) vector that includes the light quark fields, ip = (u d s)'^, and 
rh = diag (rriu, nid, rus) stands for the current quark mass matrix. For simplicity we consider 
the isospin symmetry limit, in which = = m. The currents ja^{x) are given by 

Jaix) = J d^y d'^z r{y - x) r{x - z) i){y) A„ 'ip{z) , (2) 

,rw^ /a ).%)n. A., (.). (3) 



where the form factor r{x — y) is local in momentum space, namely 

r{x-y)^l^^ e-(-^)^ r{p) , (4) 

and the matrices A^, with a — 0, ..,8, are the usual eight Gell-Mann 3x3 matrices — 
generators of SU(3) — plus Aq = •\/2/3 Isxs- Finally, the constants Tabc in the t'Hooft term 
accounting for flavor-mixing are defined by 

Tabc = ^ ^ijk ^mnl i^a)im i^bjjn i^c)kl ■ (5) 

The coupling to the Polyakov loop can be implemented by assuming that the quarks 
move in a background color gauge field (f) — iAg — ig 5^0 G^A"/2, where are the SU(3) 
color gauge fields. Then the traced Polyakov loop, which is taken as order parameter of 
confinement, is given by $ = |Tr exp(i0). In what follows we work in the so-called Polyakov 
gauge, in which the matrix is given a diagonal representation (j) = 0^A^ + 0^A®, which leaves 
only two independent variables, 0^ and 0^. 

To proceed we consider the grand canonical thermodynamical potential of the model 
within the mean field approximation. Using the standard Matsubara formalism we get 

d^p 

f,c n=—oo 



+ (6) 



. f 
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where / = {u,d,s), c = {r,g,b), and we have used the definition cu^^ = (cun - 
Un = (2n + l)nT being the usual Matsubara frequencies. The quantities 0c are defined by 
the relation = diag(0r, 0g, 4>b)- The constituent masses Sj-c are here momentum-dependent 
quantities, given by 



[uj: 



Within the stationary phase approximation, the mean field values of the auxi 



turn out to be related with the mean field values of the scalar fields cx/ by [22 



(Tu -\- G Su -\ — SuSg 



, a, + GSs + ^Sf = 



(7) 

iary fields Sf 
(8) 



The effective potential U{^,T), which accounts for the Polyakov loop dynamics, can be 
fitted taking into account group theory constraints together with lattice results, from which 
the temperature dependence can be estimated. Following Ref. [l9] we take 

1 



W($,T) 



a(T) + b{T) ln(l - 6 $^ + 8 <l>^ - 3 



rp4 



(9) 



with the corresponding definitions of a{T) and b{T). Owing to the charge conjugation 



properties of the QCD Lagrangian 



23j, the mean field value of the Polyakov loo p fi eld $ 



19|, this 



is expected to be a real quantity. Assuming that 0^ and 0^ are real-valued fields 
implies 0^ = 0, $ = [2cos(0Vr) + l]/3. 

For finite current quark masses the quark contribution to f2MFA(^) turns out to be diver- 
gent. To regularize it we follow the same prescription as in previous works j^. Namely, we 
subtract from f2MFA(^) the quark contribution in the absence of fermion interactions, and 
then we add it in a regularized form, i.e. after the subtraction of an infinite, T-independent 
contribution. From the minimization of this regularized thermodynamical potential, it is 
possible now to obtain a set of three coupled "gap" equations that determine the mean field 
values 0"^, and 0^ at a given temperature. 

We are also interested in the estimation of chiral condensates, which are given by the 
vacuum expectation values (uu) = (dd) and (ss). As usual, they can be obtained by varying 
^MFA with respect to the corresponding current quark masses. The explicit regularized 
expression for a quark condensate (//) reads 



(//) = 



d^p 
(2^ 



irif 



+ 1^ {n% + Uf,, 



(10) 



where n^^ = {1 + exp[{Ef ±i0c)/r]}^\ with Ej = yp"^ + m^. 

In what follows we analyze the chiral phase transition for a definite form factor, tak- 
ing into account the temperature dependence of effective masses, chiral condensates and 
susceptibilities. For simplicity, we have considered a Gaussian form factor 

r(/) = exp (-pV2A2) , (11) 

where A is a free parameter of the model, playing the role of an ultraviolet cut-off momentum 
scale. This parameter, as well as quark current masses and couplings in Eq. ([1]), can be 
chosen so as to reproduce the empirical values of meson properties at T = 0. We take into 



account the analysis in Ref. [22], where the parameters are fixed to obtain the empirical 
values of meson masses m^r, mx and m^/, together with the pion decay constant As 
shown in that work, in this way one can get a good description of the light pseudoscalar 
meson phenomenology at zero temperature. For definiteness we will work here with the 



parameter set GI in Ref. 221]. Namely, we use 



m = 8.5 MeV , = 223 MeV , A = 709 MeV , GA^ = 10.99 , HA^ = -295.3 . 

(12) 

The corresponding numerical results are given in Fig. 1, where we show the behavior of 
various physical magnitudes as functions of the temperature within the described nonlocal 
SU(3) model. Left (right) panels correspond to the model without (with) the inclusion of 
the Polyakov loop. In the upper panels we show the behavior of the mean field values cx^ 
and d"s (solid and dotted lines, respectively). For the sake of comparison, the curves have 
been normalized to the respective values at T = 0, which are found to be = 304 MeV and 
(Tg = 427 MeV (notice that at T = there is no effect of the Polyakov loop). It is clearly 
seen that in both cases the SU(2) chiral restoration proceeds as a smooth crossover, whereas 
there is an enhancement of the order of 80 MeV in the corresponding critical temperature 
when the effect of the Polyakov loop is taken into account. In order to properly define the 
values of these critical temperatures, we consider the chiral susceptibility Xg, defined as 

X<? = Hr— • (13) 



The phase transition temperature can be defined as the point where the susceptibility shows a 
peak, the sharpness of this peak serving as a measure of the steepness of the crossover . The 
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FIG. 1: Behavior of various physical magnitudes as functions of the temperature. Left (right) 
panels correspond to the model in the absence (presence) of a coupling to the Polyakov loop. The 
upper panels display the behavior of the mean field values aq (normalized to their value at T = 0) 
and the Polyakov loop, while central and bottom panels show the behavior of susceptibilities Xq 
and chiral condensates, respectively. 

curves showing the behavior of the susceptibilities in the nonlocal SU(3) models are plotted 
in the central panels of Fig. 1. As already mentioned, it is seen that the inclusion of the 
Polyakov loop leads to an enhancement in the critical SU(2) chiral restoration temperature, 
which brings it to a better agreement with the most recent lattice values Tcr*"^"*^ = 160 — 200 



MeV 2l|]. Our results are (see solid lines in the central panels of Fig. 1) 



Nonlocal SU(3) model without Polyakov loop : = 114 MeV 

Nonlocal SU(3) model including Polyakov loop : = 198 MeV 

In addition, it is seen that after the inclusion of the Polyakov loop there is an enhancement 
in the sharpness of the peak. The restoration of the SU(2) chiral symmetry can be also 
observed when looking to the behavior of the chiral condensates with the temperature, 
which is shown in the lower panels of Fig. 1. At T = we find (uu)^^'^ = —211 MeV 
and {ss)^^'^ = —186 MeV. As usual, the restoration of the SU(3) chiral symmetry is much 
less pronounced, due to the larger current mass of the strange quark. In fact, both in 
absence or presence of the Polyakov loop a broad peak is observed in the corresponding chiral 
susceptibility Xs- This is shown by the dotted lines in the central panels of Fig. 1, where the 
values of Xs have been normalized by a factor 20 to be comparable with the corresponding 
values of Xu- The peaks in Xs are found to be placed at about 190 (230) MeV for the model 
without (with) the Polyakov loop. It is interesting to notice that due to flavor mixing effects 
the light quark condensates (uu) exhibit a second drop at those temperatures. In addition, 
owing to flavor mixing some signal of the SU(2) chiral restoration can be noticed in the 
behavior of Xs'- one can observe a shoulder at about 115 MeV for the model without Polyakov 
loop (barely noticeable in the flgure), and a small peak at about 200 MeV for the model with 
Polyakov loop. Finally, it is also important to analyze the behavior of the mean fleld value 
of the Polyakov loop fleld $. As stated, the latter can be interpreted as an order parameter 
of the deconflnement transition, $ = and $ = 1 corresponding to conflned and deconflned 
quarks, respectively. In the nonlocal SU(3) theory under consideration, the behavior of $ 
with the temperature is shown by the dashed curve in the upper right panel of Fig. 1. By 
comparison with the behavior of the mean- fleld value d"„, it can be seen that even if the 
deconflnement transition goes smoother than the SU(2) chiral restoration, both transitions 
are found to occur at the same temperature region, located at about T = 200 MeV. 

The analysis presented above corresponds to the speciflc case of a covariant nonlocal 



e four-fermion interaction. Such an 
where it arises from the nonlocal 



c 1 
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SU(3) model, in which we have considered a separab 
interaction is motivated by the instanton liquid model 
structure of the QCD vacuum. For simplicity we have used here a gaussian-like form factor, 
choosing a deflnite parameter set to get a good agreement with meson phenomenology at 



zero temperature. In order to determine the range of validity of the conclusions extracted 
from our analysis, we have also considered different parameter sets (within the range allowed 
by meson phenomenology), and different shapes for the form factors such as e.g. Lorentzian 
functions. We have also analyzed an alternative way of including the nonlocality based on 
an effective one-gluon exchange picture of strong interactions (see Ref. 2J] for a detailed 
comparison of model features in the SU(2) case at T = 0). It comes out that the results are 
qualitatively similar in all these cases. 



To summarize, we have analyzed the role played by the Polyakov loop in the description 
of the chiral phase transition within the framework of nonlocal SU(3) chiral models with 
flavor mixing. We have found that its presence provides a substantial enhancement of the 
predicted critical temperature, bringing it to a better agreement with the most recent results 
of lattice calculations. Another interesting effect of the coupling to the Polyakov loop is that 
the phase transition becomes steeper, showing a sharper peak in the chiral susceptibility. 
We notice that, as it was done for the two-flavor case in Ref. Q], it would be of great 
interest to extend the present SU(3) model analysis beyond the mean fleld approximation, 
considering the effect of mesonic correlations. Other tasks to be addressed are the analysis 
of the behavior of meson properties with the temperature, and the extension of all these 
studies in the presence of flnite chemical potentials. We hope to report on these issues in 
forthcoming publications. 
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